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Abstract 



We introduce the commutators of lepton mass matrices to describe the phe- 



0^ ■ nomenon of lepton flavor mixing, and establish their relations to the effective 



Hamiltonians responsible for the propagation of neutrinos. The determinants 



^^ I of those commutators are invariant under matter effects, leading to an instruc- 

^^ I five relationship between the universal CP-violating parameters in vacuum 

■^ ' and in matter. In the scenario of low-energy (100 MeV < E < 1 GeV) and 

Q,' medium-baseline (100 km < L < 400 km) neutrino experiments, we illustrate 

JL , the features of lepton flavor mixing and CP violation. The terrestrial matter 

D ! effects on CP- and T- violating asymmetries in fg <-> f^ and Vg <-^ 17^ neutrino 

'"^ ' oscillations are also discussed. We demonstrate that a relatively pure signal 

of leptonic CP violation at the percent level can be established from such 

^ ■ medium-baseline experiments with low-energy neutrino beams. 
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1 Robust evidence for the long-standing anomalies of solar and atmospheric neutrinos 
has recently been reported by the Super-Kamiokande Collaboration [Q]. It strongly implies 
that neutrinos are massive and lepton flavors are mixed. The admixture of three different 
lepton flavors generally involves non-trivial complex phases, leading to the phenomenon 
of CP violation 0. Leptonic CP violation can manifest itself in neutrino oscillations. In 
reality, to measure CP- violating asymmetries needs a new generation of accelerator neutrino 
experiments with very long baselines 0]. In such long-baseline experiments the earth-induced 
matter effects, which are likely to deform the neutrino oscillation patterns in vacuum and 
to fake the genuine CP- violating signals, must be taken into account. 

Recently some attention has been paid to an interesting possibility to measure lepton 
flavor mixing and CP violation in the medium-baseline neutrino experiments with low- 
energy beams [§|-|§]. The essential idea is on the one hand to minimize the terrestrial matter 
effects, which are more signiflcant in the long-baseline neutrino experiments, and on the 
other hand to obtain the fast knowledge of neutrino mixing and CP violation long before a 
neutrino factory based on the muon storage ring is really built. Although such low-energy 
and medium-baseline neutrino experiments may somehow suffer from the problems such as 
smaller detection crosss sections and larger beam opening angles, they can well be realized 
by choosing the optimum baseline length and beam energy. They are also expected to be 
complementary to the high-energy and long-baseline neutrino experiments, towards a full 
determination of lepton flavor mixing and CP-violating parameters. 

The purpose of this paper is three-fold. First, the commutators of lepton mass matrices 
are introduced to describe the phenomenon of lepton flavor mixing, and their relations to 
the effective Hamiltonian responsible for the propagation of Dirac and Majorana neutrinos 
are established. An elegant relationship between the universal CP-violating parameters in 
matter and in vacuum, no matter whether neutrinos are Dirac or Majorana particles, can 
then be derived from the determinants of those commutators, which are invariant under 
matter effects. Secondly, we study the features of lepton flavor mixing and CP violation in 
the scenario of low-energy (100 MeV < -E" < 1 GeV) and medium-baseline (100 km < L < 400 
km) neutrino experiments. The terrestrial matter effects on the elements of the flavor mixing 
matrix and the rephasing-invariant measure of CP violation are in particular illustrated. 
Finally we analyze the CP- and T-violating asymmetries in z/g ^-> z/^ and Vg <-> z7^ neutrino 
oscillations based on a low-energy and medium-baseline experiment . It is demonstrated that 
a relatively pure signal of leptonic CP violation at the percent level can be established from 
such accelerator neutrino experiments. 

2 Let us denote the mass matrices of charged leptons and neutrinos in vacuum to 
be Ml and M^, respectively. The phenomenon of lepton flavor mixing arises from the 
mismatch between the diagonalization of Mi and that of M^, in an arbitrary flavor ba- 
sis. Without loss of generality, one can choose to identify the flavor eigenstates of charged 
leptons with their mass eigenstates. In this speciflc basis, where M/ takes the diagonal form 
Ml = Diagjme, ■ni^, r/i,-}, the corresponding lepton flavor mixing matrix V links the neu- 
trino mass eigenstates (z/i, z/2, z/3) to the neutrino flavor eigenstates (z/g, z/^, z/,-). The effective 
Hamiltonian responsible for the propagation of neutrinos in vacuum can be written as 

^eff = ^ (VMy) , (1) 



where My = Diag{mi, ■m2, WI3} with nii being the neutrino mass eigenvalues, and E ^ mi 
denotes the neutrino beam energy. When neutrinos travel through a normal material 
medium (e.g., the earth), which consists of electrons but of no muons or taus, they en- 
counter both charged- and neutral-current interactions with electrons. The neutral-current 
interactions are universal for z/e, v^ and v^. neutrinos, therefore they lead only to an overall 
and unobservable phase for neutrino mixing. The charged-current interactions are likely 
to modify the features of neutrino mixing in vacuum and must be taken into account in 
all realistic neutrino oscillation experiments. Let us use Mj, and V to denote the effective 
neutrino mass matrix and the effective flavor mixing matrix in matter. Then the effective 
Hamiltonian responsible for the propagation of neutrinos in matter can be written as 0] 
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in which Mj, = Diagjmi, 012, 1113} with rrij being the effective neutrino mass eigenvalues in 
matter. The deviation of TV^ from TicS, denoted as A?-^eff, is given by 



AK 



eff 




(3) 



in which A = \/2 G-pN^. measures the charged-current contribution to the z/gC" forward 
scattering, and Ne is the background density of electrons. Assuming a constant earth density 
profile (i.e., N^ = constant), which is quite close to the reality for most of the long- and 
medium-baseline neutrino experiments 0, one can derive the analytical relationship between 
nii and rrii as well as that between V and V from Eqs. (1) and (2) [Q. Instead of repeating 
such calculations, we hst the relevant results for rrij and V in Appendix A. Subsequently 
we concentrate on the commutators of lepton mass matrices and explore how leptonic CP 
violation in matter is related to that in vacuum. Our discussion in this section is essentially 
the extension of that in Ref. |^. 

An instructive measure of the lepton flavor mixing, i.e., the mismatch between the di- 
agonalization of Mi and that of M^ (or Mj^), is the commutators of lepton mass matrices, 
which can be defined as 
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for neutrinos propagating through vacuum; or as 
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(5) 



for neutrinos interacting with matter. In the flavor basis where Mi = Mi holds, we obtain 
X = Y, X = Y and X = Y, X = Y. It is important to note that 



MMl = VmIv^ 
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(6) 



hold in the chosen basis, no matter whether neutrinos are Dirac or Majorana particles (see 
Appendix B for a brief proof). In contrast, there is no direct relationship between M^M^, 
(or MJ,Mj,) and Hes (or ?^^). We are therefore interested only in the commutators X and 
X, which can be rewritten as 



X 
X 



VmIv^ , Ml 



(7) 



The determinants of X and X are related to leptonic CP violation in vacuum and in matter, 
respectively. To see this point clearly, we carry out a straightforward calculation and arrive 
at 



Det(X) = 2jnK-^^)n( 
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where the Greek indices run over (e, /i,r); the Latin indices run over (1,2,3); J and J are 
the universal CP- violating parameters of V and V, respectively. In vacuum J is defined 



through [^ 
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(9) 



where (a,/3,7) and {i,j,k) run over {e,fi,T) and (1,2,3), respectively. Similarly J can be 
defined in terms of the matrix elements of V. We see that the quantities Det(X) and Det(X) 
measure leptonic CP violation rephasing-invariantly. 

We proceed to find out the relationship between the universal CP-violating parameters 
J and J. In view of Eq. (6) as well as Eq. (3), we immediately realize that 
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(10) 



This interesting result implies that the commutators of lepton mass matrices are invariant 
under matter effects. As a consequence, Det(X) = Det(X), leading to an elegant relationship 
between J and J HI: 



i<j 
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(11) 



Kj 



J depends on the matter effect (i.e., the parameter A) through m^. Of course J = J if 
A = 0, and J = if J = 0. 



The results obtained above are valid for neutrinos propagating in vacuum and in matter. 
As for antineutrinos, the corresponding formulas can straightforwardly be written out from 
Eqs. (1) - (11) through the replacements V =^ V* and A =^ (—A). 

O Now we illustrate the features of lepton flavor mixing and CP violation in the 
scenario of low-energy (100 MeV < E < 1 GeV) and medium-baseline (100 km < L < 
400 km) neutrino experiments. As the large-angle MSW solution to the solar neutrino 
problem seems to be favored by the latest Super-Kamiokande data [^, we typically take 
the mass-squared difference Arn^^^ = 5 ■ 10~^ eV^. We take Aml^^ = 3 ■ 10"'^ eV^ in 
view of the present data on atmospheric neutrino oscillations. The corresponding pattern of 
lepton flavor mixing is expected to be nearly bi- maximal; i.e., {Vesl ^ 1 as supported by the 
CHOOZ experiment pf, \Vei\ ~ \Ve2\ ~ C(l), and |V^3| ~ iKsl ~ 0{l). The approximate 
decoupling of atmospheric and solar neutrino oscillations implies that one can set 
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Am2i = m2 - m^ ^ ±A"^sun , 



Of course lS.vnh,2 ~ ^^^31 holds in this approximation. Note that the phenomena of flavor 
mixing and CP violation in neutrino oscillations can fully be described in terms of four 
independent parameters of V . To satisfy the nearly bi-maximal neutrino mixing with large 
CP violation, we typically choose |Ki| = 0.816, IK2I = 0.571, |V;,3| = 0.640, and J = ±0.020 
as the four independent input parameters in the numerical calculations Q In addition, 
the dependence of the terrestrial matter effect on the neutrino beam energy is given as 
A = 2EA = 2.28 ■ 10"^ eV^E/[GeV] [|. With the help of Eq. (11) and those listed in 
Appendix A, we are then able to compute the ratios Amf-^/ Am"^^ (for i = 2,3), |Vai|/|V^i| 
(for a = e, ^,T and i = 1, 2, 3), and 3/ J as functions of E. The relevant results are shown 
respectively in Figs. 1, 2 and 3, in which Am2i > and Am'^^ > have been assumed. 
Some comments are in order. 

(a) Fig.l shows that the earth-induced matter effect on Am2i is significant, but that 
on Amli is negligibly small in the chosen range of E (i.e., Am|^ ^ ^"^31 is an acceptable 
approximation). In addition, the mass-squared differences of neutrinos are relatively more 
sensitive to the matter effect than those of antineutrinos. 

(b) Fig. 2 shows that the matter effects on |V^3| and |V^3| are negligible in the low-energy 
neutrino experiment. The smallest matrix element \Ve3\ is weakly sensitive to the matter 
effect; e.g., |Ve3|/|Ve3| deviates about 7% from unity ioi E ^ 1 GeV. In contrast, the other 
six matrix elements of V are significantly modified by the terrestrial matter effects. For 
|Ki|, |V^2| and IK-2I; the relevant matter effects of neutrinos are more important than those 
of antineutrinos. For \Ve2\, |V^i| and |K-i|, the matter effects of neutrinos and antineutrinos 
are essentially comparable in magnitude. 



^This specific choice corresponds to 612 ~ 35°, 623 ~ 40°, ^13 ss 5°, and 5 ^ ±90° in the 
standard parametrization of V ||l^ , in which J = sin 612 cos 612 sin 623 cos 623 sin ^13 cos^ ^13 sin 6 
holds. J > and J < are associated with the case of neutrino mixing and that of antineutrino 
mixing, respectively. 



(c) Fig. 3 shows that the magnitude of J may rapidly decrease with the neutrino beam 
energy E. This feature imphes that the low-energy and medium-baseline neutrino experi- 
ments are likely to provide a good chance for measurements of leptonic CP-violating and 
T-violating asymmetries. Note also that J undergoes a resonance around E ~ 100 MeV 
because of the terrestrial matter effect. 

So far we have chosen Arrigi > and Arrig^ > in the numerical calculations. As the 
spectrum of neutrino masses is unknown, it is possible that Am2i < and (or) Arrig]^ < 0. 
After a careful analysis of the dependence of J on positive and negative values of Am2i and 
Arrag^, we find the following exact relations: 
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(13) 

The validity of these equations, which are independent of both the neutrino beam energy 
and the baseline length, can easily be proved with the help of Eqs. (11), (Al) and (A2). 
The key point is that Amli and Arngj^ keep unchanged in proper arrangements of the signs 
for Ar/igi? Arn^i and A. Furthermore, there exist four approximate relations: 

J(+Am2i, +Am3i, +A) ^ J(+Am2i, -Aml^, +A) , 
J(+Am^i, +Am^i, ~A) ^ J{+Aml^, -Am^^, -A) , 
3{-Aml,, +Aml„ -A) ^ 3{-Aml,, -Amj,, -A) , 
J(-Am2i, -Amgi, +A) ^ J(-Am2i, +Aml^, +A) , (14) 

which hold to a good degree of accuracy (with the relative errors < 15%) for the chosen 
neutrino beam energy. This result implies that changing the sign of Amf^ does not affect 
the magnitude of J significantly 0. We expect that those interesting relations in Eqs. (13) 
and (14) can experimentally be tested in the near future. 

In a similar way we have carefully analyzed the behaviors of \\ai\ (for a = e, ij,,t and 
i = 1,2,3) with respect to the negative values of Am^i and Aml^. Instead of presenting 
the details of our quantitative results, we only make two qualitative remarks: (1) while J 
is not sensitive to the sign of Am|^, |Vq,j| may dramatically be suppressed or enhanced by 
changing the sign of Am"^^; (2) in contrast, the dependence of \V^ai\ on the sign of Am2i is 
less significant. 

Finally it is worth mentioning that dramatic changes of the results shown in Figs. 1-3 
do not happen, even if one allows every parameter to change in a reasonable region around 
the originally chosen value (e.g., Aml^ ~ (10^^ — 10""^) eV^, Aml^ ~ (1 — 6) • 10~^ eV^, 
iKil ~ 0.7-0.9, IK2I ~ 0.5 -0.7, 1 1/^3 1 ~ 0.5 -0.8, and \J\ ~ 0.01-0.03). This observation 



^Note that the sign of Am,32 (= ^rn'^i — A?7i2i) changes simultaneously with that of Am'^i, simply 
because |Am32| ~ |A?Ti3]^| ^ |Am,2i| holds. The sensitivity of J/ J to the signs of Arn-21 and Arn-3]^ 
has been examined in Ref. 



implies that the matter-corrected quantities of lepton flavor mixing and CP violation have 
quite stable behaviors, as a qualitative consequence of the nearly bi-maximal mixing pattern 
of V and the rough hierarchy Am'^^^ < A < Aml^j^ in the experimental scenario under 
consideration. 

4 Let us turn to CP- and T-violating asymmetries in the low-energy and medium- 
baseline neutrino oscillation experiments. The conversion probability of a neutrino z/^ to 
another neutrino z/^ in vacuum is given by [^ 

^K^^/j) = -'^T.iMVarVpjV:^V*,) ■ sin^ F,,] + 8JY[smFij, (15) 

i<j i<j 

where (a,/3) run over (e,/i), (yU, r) or (r, e), and Fij = 1.27 Amf^L/ E with L being the 
baseline length (in unit of km) and E being the neutrino beam energy (in unit of GeV). 
The probabilities of v^ -^ //„ and u^ -^ V^^ transitions can straightforwardly be read off 
from Eq. (15) with the replacement J => —J. Clearly P(I7q, -^ Va) = P^Ua -^ Vq) holds 
as a straightforward consequence of CPT invariance. The CP-violating asymmetry between 
P{i'a — > i^p) and PiVa -^ T^p) amounts to the T-violating asymmetry between Pi^a -^ ^p) 
and P(z/^ ^ z/„) 0,|T5[: 

AP = P(z/„ ^ z/^) - P{V^^V^) 
= P{v^ -^ z/^) - P{vi^ -> z/„) 
= — 16JsinP2i ■ sinPai ■ sinP32 . (16) 

Because of lArrig^l k. \Am\2\ ^ I ^^^21 1 5 ^ favorable signal of CP or T violation can be ob- 
tained only when the condition |Am2i| ~ E/L is satisfied (i.e., | sinP2i| ~ C^(l) is acquired). 
Therefore only the large-angle MSW solution to the solar neutrino problem is possible to 
meet such a prerequisite for the observation of CP violation in realistic neutrino oscillation 
experiments. 

When the terrestrial matter effects are taken into account, the probability of the Va ~-* ^a 
transition becomes 

P(z/„ -. vp) = -4 5][Re(V,,V^, V:^. V;j ■ sin^ F,,] + 8 J J] sin F,,- , (17) 

i<j i<j 

in which Fj^ = 1.27AinfjL/E. The probability P(z/^ -^ z/q,) can directly be read off from 
Eq. (17) with the replacement J =^ —J. To obtain the probability P(z7a -^ z7^), however, 
both the replacements J => —J and A => —A need be made for Eq. (17). In this case, 
P(^a "*■ ^/j) is not equal to P(z/^ — > z/^). Such a false signal of CPT violation measures the 
matter effect! Similar to Eq. (16), the CP- and T-violating asymmetries in the presence of 
matter effects can be defined respectively as 

AP = P(z/, ^ z/^) - P(z7„ ^ z7^) , 

AP = P(z/„ ^ z/^) - P{up^ v^) , (18) 

where the flavor indices (a, /3) run over (e, /x), (/x, r) or (r, e). In general, AP 7^ AP because 
of the matter-induced corrections to AP. Note that the overall matter effects residing in 



AP is expected to be tiny. Indeed AP ^ AP has numerically been confirmed to be an 
excellent approximation f^,^. The reason is simply that the matter-induced corrections 
to P(z/q, —>■ Ug) and -P(i^^ —>■ Va)-, which depend on the same parameter A, may essentially 
cancel each other in the asymmetry AP. In contrast, P(z/q -^ Va) and "PiVa -^ 1^/3) are 
associated separately with (+A) and (—A), thus there is no large cancellation of matter 
effects in the asymmetry AP. 

To illustrate, we perform a numerical analysis of the CP-violating asymmetry AP be- 
tween Ue — > i^/j, and z7e -^ T^fi transitions as well as the T-violating asymmetry AP between 
^e ~^ ^11 and v^ -^ z/g transitions. The values of the relevant input parameters are taken 
the same as before. We first choose Am\^ > and Am'^^ > to calculate AP and AP as 
functions of the beam energy E with respect to L = 100 km, 200 km, 300 km and 400 km. 
The results are shown in Fig. 4. For the chosen range of E, one can see that the matter 
effect on AP is insignificant; i.e., AP ^ AP is a good approximation. The largest deviation 
of AP from AP, of the magnitude 0.12% or so, takes place when E and L satisfy the rough 
condition L/E ~ 500 km/GeV. This observation is of course dependent on the values of the 
input parameters. We confirm that AP ^ AP holds to a better degree of accuracy than 
AP ^ AP. Following a similar analysis, we find that the CP-violating asymmetry between 
^'/i -^ ^e and u^ -^ z7g transitions amounts essentially to (— AP). It becomes clear that a 
relatively pure CP-violating asymmetry at the percent level can be obtained from such a 
low-energy and medium-baseline neutrino experiment. 

The continuous dependence of the CP-violating asymmetry AP on the baseline length 
L is shown in Fig. 5 for five different values of the neutrino beam energy E (i.e., E = 100 
MeV, 200 MeV, 300 MeV, 400 MeV and 500 MeV). Fixing E = 100 MeV, for example, 
we observe that the maximal magnitude of AP may reach 4% at L ~ 200 km and 6% at 
L ~ 280 km. Note that L should not be too large in realistic experiments, in order to keep 
the neutrino beam opening angles as small as possible [§]. For a specific experiment, one 
can select the optimum values of E and L after taking its luminosity and other technical 
details into account. 

Finally let us examine the sensitivity of AP to the signs of Am\^ and Amg^. To be 
explicit, we take L = 100 km and 400 km. The dependence of AP on E is then studied 
with respect to four possible combinations for the signs of Am^i and Amli. The numerical 
results are depicted in Fig. 6. We find that the following relations hold to a good degree of 
accuracy: 

AP(+Am2i,+Am3i) ^ -AP{-A'ml^,-Aml^) , 

AP{+Aml^, -Amgi) ^ -AP(-Am2i, +Aml^) . (19) 

To understand the above equation, we notice from Eq. (16) that AP(+Am2i, iArrig]^) and 
AP{—Amli, =i=Am3]^) have the same magnitude but the opposite signs. A slight deviation of 
AP(±A?7i2i, ±Am|^) from AP(±Am2i, ±Am|;^) arises from the terrestrial matter effects, 
which are very samll in the experimental scenario under consideration. Therefore we arrive 
at the approximate relations in Eq. (19). Note also that the sensitivity of AP to the sign 
of Am'^i is quite weak, as shown in Fig. 6. This result is certainly consistent with that for 
J in Eq. (14); i.e., changing the sign of Ar7i|^ does not affect J and AP significantly. 

5 In this paper we have introduced the commutators of lepton mass matrices to 



describe the phenomenon of lepton flavor mixing. Their relations to the effective Hamiltonian 
responsible for the propagation of neutrinos are independent of the nature of neutrinos (Dirac 
or Majorana). The determinants of those commutators are invariant under matter effects, 
leading to an elegant relationship between the universal CP-violating parameters in matter 
and in vacuum. 

We have illustrated the features of lepton flavor mixing and CP violation in the scenario of 
low-energy (100 MeV < E < 1 GeV) and medium-baseline (100 km < L < 400 km) neutrino 
experiments, In particular, the terrestrial matter effects on the elements of the lepton mixing 
matrix and on the rephasing-invariant measure of CP violation are systematically analyzed. 
Some useful relations have been found for the matter-corrected parameters of CP or T 
violation in respect to different signs of the neutrino mass-squared differences. 

We have also presented a detailed analysis of CP- and T-violating asymmetries in neu- 
trino oscillations, based on a medium-baseline experiment with low-energy neutrino beams. 
The terrestrial matter effects are demonstrated to be insigniflcant and sometimes even neg- 
ligible in such an experimental scenario. A relatively pure signal of leptonic CP violation at 
the percent level can be established from the probability asymmetry between z/^ —>■ v^ and 
z7^ — > z7g transitions, or between z/g — > z/^ and I/g — >• v^ transitions. 

To realize a low-energy and medium-baseline neutrino experiment needs high-intensity 
and narrow-band neutrino beams (e.g., 10 to 100 times more intense than the neutrino 
flux in the K2K experiment [Q). Considering conventional neutrino beams produced from 
the charged pion decay, for example, one may maximize the neutrino flux in the forward 
direction by restricting the pion beam divergence. The required radial focusing can be 
provided by a quadrupole channel and (or) magnetic horns P,p!8[]. The peak pion energy 
and energy spread within the pion decay channel are determined by the beamline optics, 
which in turn determine the neutrino spectrum. If the optics are designed to accept a large 
pion momentum spread, the resultant wide-band beam will contain a large neutrino flux with 
a broad energy spectrum; and if the optics are designed to accept a smaller pion momentum 
spread, the resultant narrow-band beam will contain a smaller neutrino flux with a narrower 
energy spectrum [|l^. A detailed study of such technical problems is certainly desirable, but 
beyond the scope of this work. 

We expect that the low-energy, medium-baseline neutrino experiments and the high- 
energy, long-baseline neutrino experiments may be complementary to each other, in the near 
future, towards a precise determination of the flavor mixing and CP-violating parameters in 
the lepton sector. 

The author is indebted to A. Romanino for her enlightening comments on this paper. 



APPENDIX A 

In the assumption of a constant earth density profile and with the help of the effective 
Hamiltonians given in in Eqs. (1) and (2), the author has recently derived the exact and 
parametrization-independent formulas for the matter-corrected neutrino masses iiij and the 
flavor mixing matrix elements V^j in Ref. [Q. The main results are briefly summarized as 
follows. 

(1) The neutrino mass eigenvalues m^ in vacuum and m^ in matter are related to each 
other through 

1 



mf = ml H — X a/x^ 



3y 



z + J3 1 - z' 



rrin 
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ml + l.-l^/^^^y 
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where x, y and z are given by [|T6l 



X = Amli + ArTig;^ + A 
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1 2x^ - 9xy + 27AAml.Aml. I Vei I 
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with Am|i and Amj^ defined in Eq. (12), and A = 2EA = 2^2 GyN^E. Of course, 
m^ = ml can be reproduced from Eq. (Al) if A = is taken. Only the mass-squared 
differences Am2i = 1112 — m^ and Anig]^ = rrig — mf are relevant to the practical neutrino 
oscillations in matter. 

At this point it is worth mentioning that m^, |V^ip, A and m^ are correlated with one 
another via an interesting equation, which has not been noticed in the literature: 



A (ml - ml) (m^j - ml) |Kfc| 
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(A3) 



where {i,j,k) run over (1,2,3) with i ^ j ^ k. Assuming |Ar7i|;^| ^ |A?7i|]^| and taking 
k = 3, one can easily reproduce the approximate analytical result for the correlation between 
m| and ml obtained in Ref. fl^ . 

(2) The analytical relationship between the elements of V in matter and those of V in 
vacuum reads: 
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in which a runs over (e, /i, r) and {i,j, k) over (1, 2, 3) with i ^ j ^ k, and 
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Obviously A = leads to \ai = Vm- This exact and compact formula shows clearly 
how the flavor mixing matrix in vacuum is corrected by the matter effects. Instructive 
analytical approximations can be made for Eq. (A4), once the spectrum of neutrino masses 
is experimentally known or theoretically predicted. 

The results listed above are valid for neutrinos interacting with matter. As for antineu- 
trinos propagating through matter, the relevant formulas can straightforwardly be obtained 
from Eqs. (Al) - (A5) through the replacements V =^ V* and A =^ —A. 



APPENDIX B 



This Appendix aims to show why M^Ml = VM^V'^ holds in the flavor basis where the 
charged lepton mass matrix is diagonal (i.e., Mi = Mi), no matter whether neutrinos are 
Dirac or Majorana particles. For simplicity, we focus on the standard charged-current weak 
interactions, in which only the left-handed leptons take part: 




- /^wcak = -^ {lye, y^, yr\ Y \ ^^ \ W+ + h.C. , (Bl) 

where the flavor eigenstates of charged leptons are identified with their mass eigenstates. As 
the lepton flavor mixing matrix V is defined to link the neutrino flavor eigenstates (z/g, ^ii-, ^t) 
to the neutrino mass eigenstates (z/i, z/2, z/3), we get 

- £„eak = ^ (^1,^2,^3)l V^l' U K + ^.C. (B2) 

in the chosen flavor basis. 

If neutrinos are Dirac particles, the mass term can be written as 




--^Dirac = {l'e,^^^l,^^T)'L Ml, U^ + h.C, (B3) 



where M^, is in general an arbitrary 3x3 matrix. It is always possible to diagonalize Mj, 
by a bi-unitary transformation: V'^MyV = M^, where V is just the flavor mixing matrix of 
Dirac neutrinos consistent with Eq. (B2). Obviously M^Ml = VM^V^ holds. 
If neutrinos are Majorana particles, the mass terms turns out to be 




- -CMajorana = " («^e, l^fi, '^t)-L M^ \ ^^t \ + h.C. , (B4) 

in which u'^ = Cv"^ with C being the charge-conjugation operator. It is well known that Mj, 
must be a symmetric matrix and can be diagonalized by a single unitary transformation: 

11 



V^M,yV* = M,y, where V is just the flavor mixing matrix of Majorana neutrinos consistent 



■2, 



with Eq. (B2) 0. Once again we arrive at M^M^ = VM^VK 

A more general neutrino mass Lagrangian involves both Dirac and Majorana terms [jl9 
In this case, one can similarly prove that the mass matrix of light (active) Majorana neutrinos 
is symmetric and satisfies M^Ml = VM^V^ . Thus Tieff = {.M^Ml) / {2E) holds in vacuum, no 
matter whether neutrinos are Dirac or Majorana particles. In matter we have an analogous 
relation between TV^^ and yijs/i\, as given in Eq. (6). 



^If the Majorana mass matrix My is diagonalized by the transformation U"^ MyU = My, one will 
see that it is U* (instead of U) linking the flavor eigenstates to the mass eigenstates of neutrinos. 
Therefore the flavor mixing matrix turns out to be ^ = U* in the chosen flavor basis. 
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FIG. 1. Ratios Am2i/A?Ti,2i and Anig^/Aml]^ changing with the beam energy E (in unit of 



GeV) for neutrinos (i^) and antineutrinos (zv), in which Ar7T,2i 
|I4i| = 0.816, and II42I = 0.571 have typicahy been input. 
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FIG. 2. Ratios |VQ,i|/|yQ,i| (for a = e,//, r and i = 1,2,3) changing with the beam energy 
E (in unit of GeV) for neutrinos {y) and antineutrinos (f), in which Attt-^i = 5 • 10~^ eV , 
Am|i = 3 • 10^3 gY2^ |y^^| ^ Q g;^g^ i^^^i ^ Q 57;^^ |y^^| ^ q g4Q^ ^^^ j ^ -to.020 have typically 

been input. 
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FIG. 4. The CP-violating asymmetries between u, 
AP in unit of 10^^) and in vacuum (b: AP in unit of 10^^) changing with the beam energy E (in 
unit of GeV), in which Am^^ = 5 • IQ-^ eV^ Am|i = 3 • lO'^ eV^ \Vei\ = 0.816, \Ve2\ = 0.571, 
|V^3| = 0.640, and J = ±0.020 have typically been input. 
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FIG. 5. The CP-violating asymmetry AP between zv, 
10~^) changing with the basehne length L (in unit of km) and the beam energy E (a: 100 McV; b: 
200 MeV; c: 300 MeV; d: 400 MeV; e: 500 MeV). Here Am^i = 5- 10"^ eV^ Am|i = S-IO"^ eV^ 
|Vei| = 0.816, \Ve2\ = 0.571, iV^sl = 0.640, and J = ±0.020 have typically been input. 
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FIG. 6. Dependence of the CP-violating asymmetry AP (in unit of 10 ^) on the beam energy 
E (in unit of GeV) and on the signs of (Am^^, Am§i) - 1: (-, +); 2: (-, -); 3: (+, +); 4: (+, 
-). Here |Amii| = 5 • 10"^ eV^ |Amii| = 3 • 10^^ eV^ \Vei\ = 0.816, \Ve2\ = 0.571, \V^3\ = 0.640, 
and J = ±0.020 have typically been input. 
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